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Abstract 



Corrections to hyperfine splitting and Lamb shift of order a^{Za)^ induced 
by the diagrams with radiative photon insertions in the electron line are 
calculated in the Pried- Yennie gauge. These contributions are as large as 
-7.725{3)a'^{Zaf/{TTn^){mr/mfm and -O.Q711{7)a'^{Za)/{7rn^)EF for the 
Lamb shift and hyperfine splitting, respectively. Phenomenological implica- 
tions of these results are discussed with special emphasis on the accuracy 
of the theoretical predictions for the Lamb shift and experimental determi- 
nation of the Rydberg constant. New more precise value of the Rydberg 
constant is obtained on the basis of the improved theory and experimental 
data. 



A steady and rapid progress in the spectroscopic measurements in recent 
years led to a dramatic increase of accuracy in the measurements of the Ryd- 
berg constant []I], 0, ground state IS* Lamb shift in hydrogen and deuterium 
[0, 0, ^ , classic 25* — 2P Lamb shift in hydrogen H, |^ , and of the muonium 
hyperfine splitting in the ground state 0, |10[ (see Table 1). 

These spectacular experimental achievements constitute a serious chal- 
lenge to the theory and intensive theoretical efforts are necessary to match 
this experimental accuracy. 

Theoretical work on the high order corrections to hyperfine splitting 
(HFS) and Lamb shift concentrated recently on calculation of nonrecoil con- 
tributions of order a^(Za)^. Their magnitude may run up to several kilo- 
hertz for HFS in the ground state of muonium, to several tens of kilohertz for 
n = 2 Lamb shift in hydrogen and may be as large as hundreds of kilohertz 
for the ground state Lamb shift in hydrogen. Contributions of such order 
of magnitude are clearly crucial for comparison of the current and pending 
experimental results with the theory. 

As was shown in for hyperfine splitting and in for the Lamb shift 
there are six gauge invariant sets of diagrams (see Fig.l), which produce 
corrections of order a'^{Za)^. All these diagrams may be obtained from 
the skeleton diagram, which contains two external photons attached to the 
electron line, with the help of different radiative insertions. All contributions 
induced by the diagrams in Figs, la — le, containing closed electron loops, 
were obtained recently in papers |ll|, [l^, |l^ for the case of hyperfine 



splitting and in papers |]T2|, |T6|, 0, for the case of the Lamb shift. 
These theoretical results are now firmly established since all these corrections 
were calculated independently by two different groups and the results of these 
calculations are in excellent agreement. 

We report below on the results of our calculation of the contributions of 
order a'^{Za)^ to HFS and Lamb shift induced by the last gauge invariant 
set of diagrams in Fig.l/. This set includes nineteen topologically different 
diagrams 0] presented in Fig. 2. The simplest way to describe these graphs 
is to realize that they were obtained from the three graphs for the two-loop 
electron self-energy by insertion of two external photons in all possible ways. 
Really, graphs 2a — 2c are obtained from the two-loop reducible electron self- 
energy diagram, graphs 2d — 2k are the result of all possible insertions of two 
external photons in the rainbow self-energy diagram, and diagrams 21 — 2s are 
connected with the overlapping two-loop self-energy graph. We have already 
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calculated contributions induced by the diagrams in Figs. 2a — 2h and Fig.2Z 



earlier Results of the calculation of the contributions produced by 

the remaining diagrams in Fig. 2 are presented below. 

Let us start with a brief description of the main features of our approach 
to calculations. As was shown in ||2^ for HFS and in [l^ for the Lamb 



shift, contributions to the energy splittings are given by the matrix elements 
of the diagrams in Fig. 2 calculated between free electron spinors with all 
external electron lines on the mass shell, projected on the respective spin 
states and multiphed by the square at the origin of the Schrodinger- Coulomb 
wave function. 

Actual calculation of the matrix elements is impeded by the ultraviolet 
and infrared divergences. Infrared problems usual more difficult to 

deal with than the ultraviolet ones. It is easy to realize that in the standard 
Feynman gauge all diagrams in Fig. 2 are infrared divergent and one has to 
introduce the radiative photon mass to regularize this divergence. Sure, the 
final result for the sum of all contributions induced by the diagrams in Fig. 2 
is infrared finite and should admit a smooth limit for the vanishing photon 
mass. However, numerical recipes used in calculations of the contributions 
to the energy shifts make it impossible to check analytically independence of 
the results on the photon mass and one has to rely on the extrapolation in 
the infrared photon mass. Of course, such approach is still feasible, but we 
have preferred to use the gauge invariance of the sum of diagrams in Fig. 2 
and to perform all calculations in the Fried- Yennie (FY) gauge for the 
radiative photons. All diagrams are infrared finite in this remarkable gauge 
and one may perform the on-mass-shell renormalization without introduction 



of the infrared photon mass (see, e.g., p2[) avoiding thus the problem of 
extrapolation to the vanishing photon mass. Of course, infrared finiteness in 
the FY gauge is not given for free, and one has to pay special attention to 
the infrared behavior of the integrand functions and to perform cancellation 
of spurious infrared divergences with the help of integration by parts over 
the Feynman parameters prior to momentum integration. 

Calculation of the contributions to the energy splittings starts with putting 
down the universal infrared diverging skeleton integrals corresponding to the 
electron line with two external photons. Each contribution of order a'^{Za)^ 
arises from radiative insertions in the skeleton graph. Corrections to hyper- 
fine splitting and Lamb shift, produced by the diagrams in Figs.l and 2 are 
given by the expressions (see, e.g. p2|, Il2|) 
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and 

TTn-^ Vtt/ Km / Jo 

where k is the magnitude of the three-dimensional momentum of the 
external photons measured in the electron mass units, = m/(l + m/M) is 
the reduced mass of the electron-muon (or electron-proton) system and Ep 
is the Fermi energy of hyperfine splitting. The functions L{k) are connected 
with the numerator structure and spin projection of each particular graph and 
describe radiative corrections to the skeleton diagram. They are normalized 
on the skeleton numerator contributions. 

It should be mentioned that some of the diagrams under consideration 
also contain contributions of the previous order in Za. Physical nature of 
these contributions is especially transparent in the case of HFS. They corre- 
spond to anomalous magnetic moment, their true order in Za is lower than 
their apparent order and they should be subtracted from the electron factor 
prior to calculation of the contributions to HFS. Analogous situation holds 
also in the case of the Lamb shift. The only difference is that this time not 
only the Pauli formfactor but also the slope of the Dirac formfactor of the 
electron is capable to produce lower order contribution to the splitting of the 
energy levels (see, e. g. ||12|, [|l^ and [0). Technically cases of lower order 
contributions both to HFS and to the Lamb shift are quite similar. Lower 
order terms are produced by the constant terms in the low-frequency asymp- 
totic expansion of the electron factor in the case of the hyperfine splitting 
and by the terms proportional to the exchanged momentum squared in the 
low-frequency asymptotic expansion of the electron factor in the case of the 
Lamb shift. 

These lower order contributions are connected with integration over ex- 
ternal photon momenta of characteristic atomic order mZa and the approxi- 
mation based on the skeleton integrals in eq.([l|) and eq.(||) is unadequate for 
their calculation. In the skeleton integral approach these previous order con- 
tributions emerge as the infrared divergences induced by the low-frequency 
terms in the electron factors. We subtract leading low-frequency terms in the 
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low-frequency asymptotic expansions of the electron factors, when necessary, 
and thus get rid of the previous order contributions. 

The results of our calculations of the contributions to HFS and Lamb 
produced by different diagrams in the FY gauge are presented in Table 2 []. 
For the total correction of order o?{Za)^ to the HFS and the Lamb shift 
produced by all diagrams in Fig. 2 we obtain 



AEg?, = -0. 6711 (7) ^ E., (3) 



and 



,(1/) _ ^ a\Zaf ^r7v^3 

i3 



AE'''> = -7. 725 (3) ^ 7 ( — Y m. (4) 



While this work was in progress two other papers were published where 



the contributions of the diagrams in Fig. 2 to HFS ||T5[ and the Lamb shift 
were calculated. The authors of these works used a completely different 
approach to calculations, in particular, they worked in the Feynman gauge 
and, hence, all contributions of individual diagrams in these works are in- 



frared divergent. The numbers cited in [^, ^ are obtained with the help 
of extrapolation of the presumably infrared finite sum of all contributions in 
Fig. 2 to the vanishing infrared photon mass. Despite the great differences in 
the approaches used in the present work and in []T3|, ^ numerical factors in 



eq.(|^) and eq.(^ are compatible with —0.63(4) in |]T5[ and with —7.61(16) in 
^ , respectively. Our numbers are about two orders of magnitude more pre- 
cise and further improvement of accuracy may be achieved. The reason for 
this increased accuracy is the use of the FY gauge, where one can avoid the 
extrapolation in the photon mass. The price we had paid for this advantage 
is the more tiresome analytic work needed to cancel all would be infrared 
divergences before integration. 

Numerically the correction to muonium HFS in the ground state produced 
by the diagrams in Fig. 2 is equal to 

A£;jJP^ = -0. 3701 (4) kHz. (5) 
and the total contribution of order a'^{Za)Ep is given by 



^Detailed account of our calculations will be presented in a separate publication. 
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/^E^Hp's^^ = 0- 4264 (4) kHz. (6) 

Taking into account other theoretical contributions to HFS and especially 
some small contributions obtained recently (see, e.g., reviews in |]T3|, and 
using for calculation the value of the fine structure constant as obtained in 



2^] one may obtain the theoretical value for the muonium HFS in the ground 



state 

^Ehfs = 4 463 302. 55 (0.18) (0.18) (1.33) kHz, (7) 

where the first error in parenthesis reflects the uncertainty of the fine 
structure constant itself and the second is induced by the uncertainty of the 
contribution of order alZaYEp. The third, and by far the largest contribu- 
tion to the error in the theoretical value of HFS is defined by the experimental 
error in measuring electron- muon mass ratio m/M. 

The agreement between theory and experiment is excellent. We will not 
dwell on the HFS problem any more here since the phenomenological situa- 
tion and the influence of the result in eq. (|^) on the value of the electron-muon 
mass ratio and the fine structure constant was discussed in great detail re- 
cently |T3|, 



The case of Lamb shift deserves more comments. Numerically the correc- 
tions to the IS and 25" Lamb shifts produced by the diagrams in Fig. 2 are 
equal to 

AEi^^^(lS) = -334. 2 (1) kHz, (8) 
AE'i^\2S) = -41. 78 (2) kHz, 
while respective total contributions of order a'^{Za)^m are given by 

Ae[^""^^^(1S) = -296. 9 (1) kHz, (9) 
^^aa-i/)^25) = -37. 12 (2) kHz. 

Let us discuss the sign and the scale of the correction in eq.(||). The sign 
may be determined by considering the electron factor, as defined in eq.(^). 
The low-frequency asymptotic behavior of the electron factor is described by 



the expresssion (see, e.g., |T7 
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L{k) ^ (-2Fi'(0) - ^^2(0)) = -0.7046225 k^, (10) 

where Fi{k'^) and F2{k'^) are the two- loop contributions to the ordinary 
Dirac and Pauh form factors of the electron (contribution of the graphs with 
vacuum polarization insertions in the photon line is omitted in eq.(|TU|)), re- 
spectively. We use in eq. (p!OD the well-known values for the slope of the Dirac 



form factor and of the Pauli form factor at zero p9|, |30[. As was explained 
above one has to subtract from the electron factor this leading low-frequency 
term which produces contribution to the Lamb shift of previous order in Za. 

It is well known from the general principles that the unsubtracted elec- 
tron factor has at most logarithmic behavior at infinity. Hence, the high 
momentum behavior of the subtracted electron factor is completely defined 
by the subtraction term in eq. (p!OD . Then it is clear from eq.(0), where the 
subtracted electron factor plays the role of the integrand, that the contri- 
bution to the Lamb shift induced by the graphs in Fig. 2 has the negative 
sign. One may even make an estimate of this contribution from the known 
asymptotic behavior of the integrand but we choose a less technical path in 
discussion of the magnitude of this contribution. 

It may seem at first sight that the magnitude of the corrections induced 
by the diagrams in Fig.l/ as presented in eq.(|^) are too large. We would 
like to emphasize that, quite opposite, this correction has exactly the scale 
one had to envisage before calculations. Let us discuss this point in shghtly 
more detail. It is helpful to recollect that the main contribution to the Lamb 
shift is a radiative correction itself and so it is misleading to normalize all 
contributions to the Lamb shift with the help of this leading order contri- 
bution. In this respect the case of the Lamb shift differs drastically from 
the case of HPS, where the leading Permi contribution is not the radiative 
correction but the classic effect of the interaction of two magnetic moments 
and sets the natural scale for all radiative corrections. Main contribution 
to the Lamb shift has the form 4:m{Za)'^/n^ x slope of the Dirac formfactor, 
where the slope is roughly speaking a/ir ■ (1/3) ■ ln(Za)~^. The skeleton 
factor which sets the scale for the different contributions to the Lamb shift 
is Am{Za)'^/n^ and to make an estimate of any correction to the Lamb shift 
one has to extract this skeleton factor. All other entries in the leading or- 
der contribution to the Lamb shift are produced by the radiative correction, 
and it is necessary to take into account that the number which should be of 
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order one, as predicts the common wisdom for radiative corrections, is the 
factor 1/3 before the logarithm which remains after extraction of the factor 
a/vr, characteristic for the one- loop radiative corrections. The other subtlety 
to be taken into account in estimating the orders of magnitude of different 
corrections is that, unlike the case of radiative corrections to the scattering 
amplitudes, in the bound state problem not every factor a is accompanied 
by an extra factor tt in the denominator. This is a well known feature of the 
Coulomb problem. 

Let us consider as an exercise in the art of making educated estimates the 
correction of order aiZaf'm calculated analytically long time ago |]2^, ^ . 
According to the considerations above the scale of this correction should be 
set by the factor 4a(ZQ;)^/n^m and the only problem of the theory is to 
calculate the number of order one before this factor. Analytic calculation 
[0, |2§ produces this factor in the form 1 + 11/128 - l/21og2 ^ 0.739 in 
excellent agreement with our qualitative considerations. Now it is easy to 
realize that the natural scale for the correction of order a^i^Zaf is set by the 
factor 4a^(Za)^/(7™^)m. The coefficient before this factor obtained above 
and in is about —1.9 and there is nothing unusual in its magnitude for 
a numerical factor corresponding to a radiative correction. 

Consider now current status of the Lamb shift theory. Theoretical pre- 
dictions presented below are obtained with the help of the expressions for 



the Lamb shift contributions as collected in the reviews [31, 32], amended. 



besides corrections obtained above and in with some other recent results 
presented in the Table 3. Note that the correction of order oP'iZaf' in this 
Table is again of reasonable magnitude since its scale is set by the factor 
4a^(Za)^/(7r^n^)m as one may easily check with the help of the arguments 
used above in the discussion of the contribution of order c^i^Zaf . On the 
background of this factor numerical factor 2/27 before the logarithm cube is 
quite moderate. 

Last line in Table 3 contains a new recoil correction corresponding to the 
insertions in the Coulomb photon of the muon or hadron vacuum polariza- 
tion operator^. Respective contribution for the muon insertion contains an 
evident extra electron-muon mass ratio squared suppression factor relative 
to the leading vacuum polarization contribution. We estimated hadron con- 

^K. Pachucki and S. Karshenboim are also considering these contributions (private 
communication from S. Karhenboim). 
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tribution approximating the spectral function below 1 GeV according to the 
vector dominance model and above 1 GeV we simply used the asymptotic 
quark value for the spectral function^. 

We use in calculation of the theoretical values for the Lamb shift new 
values for the self-energy contributions to the coefficient Cqq for the 15*- 
and 25'-states and to the function Gse for the 2Pi/2-state and for 



the 4S'i/2-state ^8[. We also use new values |^ Gvpi^Si/2) = —0.6187, 
Gvpi2Si/2) = -0.8089, G'yp(2Pi/2) = -0.0640 and Gvp{4:Si/2) = -0.8066 
for the Uehling part of the vacuum polarization contribution. These numbers 
are the sum of contributions of order a{Za)^ and of additional terms of higher 
order in Za. 

From the theoretical point of view the accuracy of calculations is limited 
by the magnitude of the yet uncalculated contributions to the Lamb shift. 
First, there are pure recoil contributions of order {Za)^{m/M)m which may 
be as large as 1.3 kHz for the 15'-state and 0.16 kHz for the 2S'-state. As far as 



we know, such terms were never discussed in the literature Il31| . It is claimed 
in 



m 



that there are no such contributions besides recoil corrections obtained 
0, but, unfortunately, no proof of this statement was published. We will 
accept that such corrections are absent, especially taking into account that 
according to the estimates above they are in any case too small to influence 
really the comparison of the theory with the experimental results. 

Unknown correction of order a^{Za)^ is induced by the three-loop slope 
of the Dirac formfactor of the electron and by the three-loop electron vac- 
uum polarization. Natural scale for this correction is set by the factor 
4:a^{Za)^/{7r^n^)m and we envisage the contributions about 17 kHz for the 
IS-state and 2 kHz for the 2S'-state. 

Next come uncalculated corrections of order a'^{Za)^. Contribution of 
this order is a polynomial in ln(Za)~^, starting with log cube. The factor 
before log cube was calculated in |Q and the contribution of the log squared 
terms to the difference El{1S) — 8El{2S) was obtained in However, the 
calculation of respective contributions to the separate energy levels is still 
missing. In this conditions it is fair to take the log squared contribution to the 
interval El{IS) — 8El{2S) as an estimate of the scale of all yet uncalculated 
corrections of this order. We thus assume that uncertainties induced by the 
yet uncalculated contributions of order a'^{Za)^ constitute 15 kHz and 2 kHz 

•^The derivation of this result will be published elsewhere. 
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for the IS"- and 2S'-states, respectively. 

The scale of the self-energy correction of order a{Zay is set by the factor 
4a{Zay /n^ . This contribution is a linear polynomial in ln(Za)~^. We are 
aware about two recent attempts to make an estimate of this contribution 
0, 1^, but, unfortunately, its final magnitude seems to be still unavailable^. 
Relying on the scale factor above and the estimates in |^ we assume that 
the corrections of order a{Zay are as large as 17 kHz and 2 kHz to the IS- 
and 2S'-states, respectively. 

All other theoretical contributions to the Lamb shift are smaller than 
those just discussed. Hence, we assume that the theoretical uncertainty of 
the expression for the Lamb shift is about 28 kHz for the IS'-state and about 
4 kHz for the 2S'-state. 

The other limit on the accuracy of the theoretical calculation of the Lamb 
shift is put by the accuracy of the measurements of the proton rms charge 
radius. As is well known there are two contradictory experimental results 
for this radius [^l], ^ and at least one of these experimental results should 
be in error. The accuracy of the proton rms charge radius claimed by the 



authors of |^ ^ produces uncertainty about 32 kHz for the 15'-state and 
about 4 kHz for the 2S'-state. 

Let us compare theoretical and experimental data for the classic 25*1/2 ~ 
2Pi/2 Lamb shift. The most precise experimental data as well as the results of 
our theoretical calculations are presented in Table 4. Theoretical results for 
the energy shifts in Table 4 contain errors in the parenthesis where the first 
error is determined by the yet uncalculated contributions to the Lamb shift, 
discussed above, and the second reflects the experimental uncertainty in the 
measurement of the proton rms charge radius. We have used experimental 
result taking into account recent theoretical correction discovered in [0]. 
Note, however, that this correction does not effect any of our conclusions. 
There are two immediate conclusions of the data in Table 4. First, as already 
mentioned in , the results of the proton rms radius measurement in |^ 



should be in error since respective value of the proton charge radius is clearly 
inconsistent with all results of the Lamb shift measurements. Second, we have 
to reject either the result of the most precise measurement of the 2S'i/2— 2Pi/2 
splitting, or the experimental value of the proton charge radius as measured in 



^P. Mohr is now working on the extraction of this correction from his respective high 
Za results (Private communication from P. Mohr). 
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Q since the Lamb shift value in contradicts theoretical value calculated 
employing the rms radius in [4^] by more than five standard deviations. 
Results of two other measurements of the classic Lamb shift are compatible 
with the theory, so we will below accept the value of the proton charge radius 
as obtained in |^ . We will return to the numbers in three last lines in Table 



4 below. 

We do not include theoretical predictions for the deuterium Lamb shift in 
Table 4, since, taking into account current discrepancies in determination of 
the deuteron charge radius and solid status of the Lamb shift theory, it seems 
preferable to use the deuteron Lamb shift data for extracting the value of 
this charge radius rather than for comparison of the Lamb shift theory and 
experiment. 

Next we turn to the discussion of the IS" Lamb shift. Unfortunately, its 
extraction from the experimental data is less straightforward. The experi- 
mentalists managed to separate measurement of the IS" Lamb shift from the 
measurement of the Rydberg constant by comparing the frequencies of two 
transitions with different main quantum numbers and excluding the large lev- 
els separation depending on the Rydberg constant. In this manner certain 
combinations of IS", 25" and of higher levels Lamb shifts are experimentally 
obtained. It is pretty easy to compare these experimental data in ^, ^, Q 
with the theory above and after trivial calculations we have found an excel- 
lent agreement between the theory and experiment. We will not put down 
these results here. 

Unbiased extraction of the IS" Lamb shift from the experimental data is 
still a problem. It is impossible to avoid using to this end the experimental 
value of the 25 Lamb shift, and the emerging value of the IS* Lamb shift 
depends thus on the experimental result for the 25*1/2— 2Pi/2 splitting. Higher 
levels Lamb shifts, which also enter the problem may be safely calculated 
with sufficient accuracy. The standard approach accepted by all experimental 
groups consists in adopting one or the other 25i/2 — 2Pi/2 experimental result 
and extracting thus the value of the 15 Lamb shift. All values in Table 1 for 
the 15 Lamb shift are obtained in this manner with the help of experimental 
values in or in for the classic Lamb shift. These values should be 
compared with our theoretical prediction 

AEl{1S) = 8 172 729 (28) (32) kHz, (11) 
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where again the first error is determined by the yet uncalculated contribu- 
tions to the Lamb shift and the second reflects the experimental uncertainty 
in the measurement of the proton rms charge radius. 

The results of all experiments mentioned in Table 1 are pretty consistent 
and their agreement with the theoretical value in eq.(|TTD is satisfactory but 
not spectacular. It is necessary to recollect at this point that the "experimen- 
tal values" in the Table depend on the experimental value of the 25*1/2 — 2Pi/2 
Lamb shift adopted in their extraction and the change of the IS* Lamb shift 
value under transition from one experimental 251/2 — 2Pi/2 Lamb shift value 
to another may be significant. 

It would be helpful to find a way to extract the value of the IS* Lamb 
shift from the experimental data unambiguously without reference to the 
251/2 — 2Pi/2 experimental results, which, while being compatible with the 
theory, seem to be somewhat larger than the theoretical prediction. A nat- 
ural way to reach this goal is to use the theoretical relation between the 15" 
and 25* Lamb shifts. A good deal of theoretical contributions to the Lamb 
shift scale as and, hence, vanish in the difference 8El{2S) — Ei{lS). In 
particular, all main sources of the theoretical uncertainty, namely, proton 
charge radius contribution and almost all yet uncalculated corrections to the 
Lamb shift mentioned above vanish. Significant contribution to this differ- 
ence may be produced only by the term of the form a'^{Za)^ \n{Za)~'^, which 
was calculated recently . This term produces correction about 14 kHz and 
the accuracy of the difference under consideration is determined by the yet 
uncalculated single log contribution of the same order. Such term would 
not change the log squared term by more than fifty percent and, hence, the 
uncertainty of the difference under consideration is about 7 kHz. Hence, we 
obtain the relation 



8El{2S) - El{IS) = A, (12) 
where A = 187 234 (7) kHz. 

Now one may obtain self-consistent values for the IS Lamb shift directly 
from the experimental data in Refs.[^ ^ with the help of the relations 

8 32 ^ 

El{1S) = -(Fi5-25 - 4F25-45,/J " yi?L(45i/2) + -A, (13) 

8 32 5 

El{1S) = -{F,s-2S - ^F2s-AP,/,) - —El{4Pi/2) + -A, (14) 
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El{1S) = J^{F^s-2S - yi^25-8D,/J - —E^iSD^/^) + -A. (15) 

Numerical results are presented in Table 5. These results have somewhat 
larger errors than the respective results in Table 1, however, they do not 
depend on the experimental value of the 2S'i/2 — 2Pi/2 Lamb shift and on 
the value of the proton charge radius. The accuracy of the self-consistent 
numbers in Table 5 is mainly determined by the accuracy of the frequency 
measurements in ^ Factor 4-5 reduction of the experimental errors 
would lead to a self-consistent determination of the 15 Lamb shift with the 
same accuracy as for the values cited in Table 1. One may even invert 
the usual approach to the 2Si/2 — 2Pi/2 and IS' Lamb shift values and ex- 
tract values of the 2S-2P Lamb shift from the respective self-consistent IS* 
values (see three last lines in Table 4). These values of the 2S'i/2 — 2Pi/2 
Lamb shift are consistent with the results of the direct measurements of the 
25*1/2 ~ 2Pi/2 Lamb shift but have somewhat larger error bars. However, 
self-consisitent values of the 25i/2 — 2Pi/2 Lamb shift would become quite 
competitive with the results of the direct measurements after the 4-5 times 
reduction of the current experimental errors in the frequency measurements 
would be achieved. 

Reduction of the errors of the values of the 15 and 25i/2 — 2Pi/2 Lamb 
shifts opens new ways to a more precise determination of the Rydberg con- 
stant. We would like to mention two new directions in the determination of 
the Rydberg constant value besides the one adopted now (see, e.g., [l|, ||]). 
First, one can use the self-consistent value of the 15 Lamb shift and respec- 
tive 25i/2 — 2Pi/2 Lamb shift to get the value of the Rydberg constant. Today 
such approach leads to a loss of accuracy in comparison with the current ex- 
perimental value of the Rydberg constant (see Table 6, where the first error 
in the self-consistent values of the Rydberg constant is determined by the ac- 
curacy of the self-consistent Lamb shift values and the second is determined 
by the accuracy of the frequency measurement), but greater accuracy may 
be achieved in future. Important advantage of such approach is that the 
value obtained in this way is independent of the direct experimental results 
on 25i/2 — 2Pi/2 Lamb shift and of the value of the proton charge radius. 
Second new approach is simply to reject the experimental data on the Lamb 
shifts and to use for the determination of the Rydberg constant directly the 
data on the frequencies of transitions between the levels with different main 
quantum numbers. Such approach becomes feasible now since the accuracy 
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of the theoretical expression for such transitions is defined by the theoretical 
error of the expression for the IS" (or 25") Lamb shift which is about 28 kHz 
(and is even smaller for the 25" Lamb shift) as discussed above and is thus 
smaller than the experimental error of the frequency determination. Respec- 
tive values of the Rydberg constant are again presented in Table 6, where the 
first error is determined by the accuracy of the theoretical expression, the 
second is defined by the experimental error of the frequency measurement, 
and the third one is determined by the experimental error in the determina- 
tion of the proton charge radius. The values of the Rydberg constant in two 
last lines in Table 6 derived from independent experimental data 0, |1| are 
pretty consistent. These values are more accurate than the ones obtained 
by other methods and are the most precise contemporary values of this con- 
stant. Natural drawback of this approach is, of course, the dependence of 
the obtained value of the Rydberg constant on the proton charge radius. 

In conclusion we would like to emphasize that the high accuracy of the 
Lamb shift theory opens new perspectives in determination of the Rydberg 
constant and of the Lamb shift in the IS*- and 2S'-states. Four directions of 
experimental investigations, namely, more precise measurement of the tran- 
sitions between levels with different main quantum numbers, more precise 
measurement of the IS" and 25* Lamb shifts, and direct measurement of the 
proton charge radius seem especially promising. It is very important that all 
these experiments are mutually complementary, since they may lead to the 
values of the Rydberg constant of comparable accuracy based on the different 
kinds of experimental data. On the theoretical side calculation of the still 
unknown corrections to the energy levels discussed above with the goal of 
reduction of the theoretical error in determination of the IS* Lamb shift to 
the level of 1 kHz (and, respectively, of the 25* Lamb shift to several tenth 
of kHz) seems to be both quite perspective and feasible. 
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Table 1. Experimental Results 



Interval 


(kHz) 


Hydrogen, lS'1/2 [i 


8 172 860 (60) 


Hydrogen, lS'1/2 


8 172 815 (70) 


Hydrogen, lS'1/2 [i 


8 172 844 (55) 


Hydrogen, 2Si,2 - 2Pi/2 i 


1 


1057 845 (9) 


Hydrogen, 2Si/2 - 2Pi/2 i 




1057 857.6 (2.1) 


Hydrogen, 2Si,2 - 2Pi/2 i 


1 


1057 839 (12) 


Muonium, HFS |[ 


4 463 302.88 (16) 
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Table 2. Corrections to HFS 



and Lamb shift 





HFS 


Lamb shift 

irn' ^ m ' 


a 


9/4 





b 


-6.65997(1) 


2.9551(1) 


c 


3.93208(1) 


-2.2231(1) 


d 


-3.903368(79) 


-5.238023(56) 


e 


4.566710(24) 


5.056278(81) 


f 


-3.404163(22) 


-1.016145(21) 


9 


2.684706(26) 


-0.1460233(52) 


h 


33/16 


153/80 


i 


0.05524(21) 


-5.51680(87) 


3 


-7.14860(39) 


-7.76648(79) 


k 


1.465834(20) 


1.959589(33) 


I 


-1.983298(95) 


1.74815(38) 


m 


3.16956(16) 


1.87541(49) 


n 


-3.59566(14) 


-1.30626(49) 





1.80491(45) 


-12.0641(16) 


P 


3.50608(16) 


6.13527(90) 


q 


-0.80380(15) 


-7.52272(83) 


r 


1.05298(18) 


14.3622(15) 


s 


0.277203(27) 


-0.930291(78) 


Total 


-0.6711(7) 


-7.725(3) 
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Table 3. Some New Contributions to the Lamb Shift 



Level 


AE 






nSy2 34] 


(41n2-|)(^^m 


nS'i/2[ie] 


y 9 27/ M TT^nS I „ J 




rr/2iV ,3 14 , i V , 3 
1LV480 ^ 16n 15n2 On^ "^^V192 ^ 32n 
+ ^10.328 


nS 1/2 





Table 4. 25i/2 - 2Pi/2 Lamb Shift 



Source of the value 


(kHz) 


Experimental result [^] 


1 057 845 (9) 


Experimental result 


1 057 857. 6 (2.1) 


Experimental result [|] 


1 057 839 (12) 


Theory, = 0.805 (11) fm [|1|] 


1057 810 (4) (4) 


Theory, Vp = 0.862 (12) fm [0 


1057 829 (4) (4) 


Self-consistent IS |P] 


1057 854 (16) 


Self-consistent IS 0, [l|] 


1057 835 (15) 


Self-consistent IS 0] 


1057 847 (13) 



Table 5. Self-Consistent IS Lamb shift values 



Source of the value 


AE (kHz) 


Ref. 


i 


8 172 915 (129) 


Ref. 


i 


8 172 763 (117) 


Ref. i 


8 172 858 (107) 


Theory, this work 


8 172 729 (28) (32) 
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Table 6. Rydberg Constant 



Source of the value 


Roo (cm"^) 


Experiment |Q 


109 737. 315 684 1 (42) 


Experiment 


109 737. 315 683 4 (24) 


Self-consistent 


ii 


109 737. 315 686 8 (58) (20) 


Self-consistent 


EE 


109 737. 315 681 1 (52) (14) 


Theory and [0] 


109 737. 315 679 7 (12 )(20) (14) 


Theory and 




109 737. 315 680 2 (05) (14) (06) 
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Figure Captions 

Fig. 1. Six gauge invariant sets of graphs producing corrections of order 
a\Zaf. 

Fig. 2. Nineteen topologically different diagrams with two radiative 
photon insertions in the electron hne. 
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